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Crystalline semiconductors may exist in difterent polytypic phases with significantly diff'erent 
electronic and optical properties. In this paper, we calculate the electronic structure and optical 
properties of diamond. Si and Ge in the lonsdaleite (hexagonal-diamond) phase. We use an empirical 
pseudopotentials method based on transferable model potentials, including spin-orbit interactions. 
We obtain band structures, densities of states and complex dielectric functions calculated in the 
dipole approximation for light polarized perpendicular and parallel to the c-axis of the crystal. We 
find strong polarization dependent optical anisotropy. Simple analytical expressions are provided 
for the dispersion relations. We find that in the lonsdaleite phase, diamond and Si remain indirect 
gap semiconductors while Ge is transformed into a direct gap semiconductor with a significantly 
smaller band gap. 



PACS numbers: 71.15.Dx, 71.20.-b, 78.20.-e, 78.20.Ci 
I. INTRODUCTION 

It is well known that under extreme conditions a crys- 
talline material may undergo a structural transition to 
a phase which remains effectively stable upon return- 
ing to standard thermodynamic conditions. The exis- 
tence of such phases, polymorphism in compounds and 
allotropy in elemental crystals, results in materials with 
differing electrical and optical properties, such diamond 
and graphite. Polytypism is a particular case of poly- 
morphism in which the coordination number does not 
change. For example, III-V and group-IV semiconductors 
can crystallize in their non-naturally occurring polytypic 
phase while maintaining their tetrahedral coordinatiorP. 
For example, III-V semiconductors can crystalizc in ei- 
ther the zincblende (ZB) phase or the hexagonal wurtzite 
(WZ) phase. Under extreme conditions cubic diamond 
(CD) transforms to a hexagonal wurzite crystal, known 
as lonsdaleite (LD). The form of diamond was not dis- 
covered until 1967 when it was found in a meteorit^l It 
was synthesized soon thereaftei^ and is now believed to 
be the hardest known substance^ . 

The high pressure phases of Si and Ge have been 
investigatecP'^ and it has been determined that both un- 
dergo a series of structural phase transitions from cu- 
bic, to /3-Sn, to simple hexagonal, to an orthohombic 
phase, to hexagonal closed packed, to a face centered cu- 
bic phas^^. On the other hand, III-V semiconductors 
typically crystalize either in /3-Sn, nickeline (NiAs) or 
rocksalt structures when subjected to high temperature 
and pressur^i^. More recently, the growth of WZ phase 
bulk GaAs was achieved under extreme condition^. 

Extreme temperatures and pressure are not the only 
ways to achieve WZ or LD growth. Laser ablation 
have been used to synthesize stable LD phase Si^ III- 
V nanowires tend to crystalize in the WZ phas^^M^l 
This has been attributed to various factors such as the 



small nanowire radiP^'^, growth kineticJ^^l^ interface 
energieJ^and electron accumulation at the catalyst's in- 
terstitial site^^. It is now believed that the tendency of 
nanowires to crystalize in WZ/LD phase may be true 
for group-IV semiconductors as weh-'^^. It has been ex- 
perimentally found that Si nanowires with a radius in 
excess of 10 nm tend to crystalize in the LD phas^^ and 
a number of recent theoretical investigations have con- 
firmed that the hexagonal LD phase is the more stable for 
Si nanowires exceeding certain critical dimensionii^ll20| 
Similar structural phase transitions are expected for Ge 
nanowires as welP^. 

Semiconductor nanowires have attracted much inter- 
est due to their potential applications such as photo- 
voltaic cellj^^ J^, n ano-mechanical resonator array^^ll^ 
THz detectorp Ea si ngle photon detector^^aMl^ field- 
effect transistor^^Smi^ single-electron transistors, and 
other device^^^ll2ll in addition, the quasi 1-D nature of 
nanowires allows materials with large lattice mismatches 
to be combined to form hetrostructures that are not pos- 
sible in planar structures. The design and characteriza- 
tion of such devices requires an understanding of the elec- 
tronic and optical properties of WZ/LD phase semicon- 
ductors. Even though bulk LD phase Si, Ge and diamond 
have been synthesized in the laboratory, their electronic 
structure still remains experimentally unverified. 

Since these materials may be used for spin based de- 
vices, spin-orbit interactions should be included. The 
band structures and dielectric functions for LD phase 
diamond, Si and Ge has been calculated using empiri- 
cal pseudopotential^, however these calculations were 
done without the inclusion of spin-orbit interactions. 
The bandstucture of LD phase diamond Si and Ge have 
also been calculated using density functional theorjESlEI]^ 
however these methods have well known shortcomings in 
predicting energies of excited states. Therefore, there 
is a need for accurate empirically based band structures 
including spin. 
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In this paper, we first present bulk electronic band 
structure calculations for diamond, Si and Ge in the LD 
phase using empirical pseudopotentials with the inclu- 
sion of spin-orbit interactions. These calculations are 
based on transferable model pseudopotentials assuming 
an ideal LD structure. The spherically symmetric ionic 
model potentials are first obtained by fitting the calcu- 
lated bulk energies of the cubic polytype to experiment. 
The band structure of the LD polytype is then obtained 
by transferring the spherically symmetric model pseu- 
dopotentials to the LD crystal structure. These model 
potentials are expected to be transferable between poly- 
types due to the similarities in their crystal structures. 
Like the cubic structure, the LD structure is built from 
tetrahedrons of the same atom but are stacked differently. 
In both structures all of the nearest neighbors and nine 
out of the twelve second nearest neighbors are at identi- 
cal crystallographic locations'*^ and all the second nearest 
neighbors are equidistant. Hence the local electronic en- 
vironment should be very similar in both polytypes. This 
method has proven to be quite successful in obtaining the 
bulk band struc tures of various semiconductor polytypes 
in the pasPsEMsSl, have used this method to predict 
the band structures of WZ phase III-V semiconductor^^ 
and our calculations are in excellent agreement with ex- 
periment for the cases for which the WZ band gaps are 
known (namely GaAs, InP and InAs). Since then, a num- 
ber of recent experiments have provided further confor- 
mations of the predicted band gaps (and their respective 
symmetries p21l56] ^]^g effective mass of InAs^,jpl. 

We have also calculated the dielectric functions for di- 
amond. Si and Ge in the LD phase in the linear response 
regime within the electric dipole approximation for light 
polarized parallel and perpendicular to the c-axis of the 
crystal. These calculations are in part motivated by re- 
cent experimental results showing that the photolumines- 
cence(PL) intensity in nanowires is strongly polarization 
dependentP^miil^ including Si nanowires^. Nanowire 
polarization effects can arise from the dielectric mismatch 
at the wire surface'^, electron-hole recombination selec- 
tion rules®"* as well as from the underlying WZ-type 
crystal structur^^^HZD xhe dominant effect could be 
more easily identified if the optical dielectric functions 
for these semiconductors were easily available. However, 
experimentally measuring optical dielectric functions for 
LD phase semiconductors is difficult due to the extreme 
conditions required for growing bulk samples. 

Our dielectric function calculations are carried out 
within the one-particle picture. Much has been done in 
recent years to include the effects of two particle con- 
tributions (such as electron-hole interactions) in the di- 
electric functiorP^t^. In the case of first principles cal- 
culations, this improves agreement with experiment for 
the low energy part of the optical spectrurrP^l. How- 
ever, the dominant contribution to the dielectric function 
comes from one-particle calculations, which by itself tend 
to be in good agreement with experimenlPSHU]. More- 
over, while two-particle corrections improve ab-initio di- 



electric functionJ^, empirical methods include some two- 
particle effects through their fit to the experimental data 
(which includes such effects). Typically, dielectric func- 
tions calculated using EPMs are in good agreement with 
experiments*^. 

The required momentum matrix elements, required for 
our calculations, are obtained from the pseudopotential 
wave functions. In general, momentum matrix elements 
for pseudopotential calculations need to be corrected for 
the missing the core states. Such corrections are typically 
done using nonlocal term^SSMI^ which can also_account 
the exchange and correlations effects as welP^tSal, jn our 
calculations, nonlocal effects are included in the form of 
spin orbit interactions only. 

Another problem with ab-initio dielectric function cal- 
culations is that £0 is often overestimatecP^HiH This 
is then improved upon by the inclusion of local field 
effect such as electron-phonon interactions, which typ- 
ically affect the low frequency part of the dielectric func- 
tion (terahertz regime) . However, this is not necessary in 
our case since local field effects shift the peak position^^. 
The EPM includes such effects through the fitting to ex- 
periments which necessarily include the effects. 

We instead adopt a simpler approach to take the miss- 
ing core states and local field effects into account. We 
correct the static dielectric function, eg, for the unknown 
polytype by making use of the known eo of the cubic 
phases of diamond. Si and Ge. First, the optical dielec- 
tric functions for these cubic group-IV semiconductors 
are calculated. The optical sum rules are then used to 
obtain a set of constants which normalizes the calculated 
cubic eo to their respectively known experimental val- 
ues. Since the constituent element of each polytype is the 
same, corrections to account for the missing core states 
should be nearly the same and transferable between poly- 
types. These normalization constants are then used to 
correct eo for the unknown polytypes(LD), which there- 
fore also corrects the LD phase dielectric functions as 
well. Unlike calculations involving local field effects, our 
simple approach fixes eo without shifting peak positions. 

This paper is organized as follows. Section [TT| outlines 
the empirical pseudopotential method, followed by a de- 
scription of both the cubic and hexagonal band struc- 
tures, and a discussion of the respective symmetries of 
the two polytypes in section [TTll Our results are in section 
|IIIB| where we present the calculated band structures, 
densities of states (DOS), effective masses for zone center 
states and transition energies at various high symmetry 
points . The optical dielectric functions and reflectivity 
spectra are given in section |IV[ Finally, we summarize 
our results in section IVl 



II. EMPIRICAL MODEL PSEUDOPOTENTIALS 

We use the empirical pseudopotential formalism of Co- 
hen and Chelikowsky^. However, rather than discrete 
form factors we use continuous model potentials so they 
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are transferable between polytypes. The pseudopotential 
Hamiltonian consists of the kinetic, local pseudopotential 
{Vpp) and SO interaction {Vso) terms, 



H 



2m 



+ Vpp + Vso. 



(1) 



In a periodic crystal Vpp can be expanded in terms of 
plane waves as 



Vpp{ 



N 



(2) 



where G are reciprocal lattice vectors, a labels the atom 
type, V^^(G) is the form factor of the a*'' type of atom, 
N is the number of atoms per unit cell of a given type, 
and Ta.i is the position of atom number j of type a. 

For a compound with only one type of atom, the pseu- 
dopotential is simply 

(3) 



{G'\Vpp\G) 



y^^(G'-G)5(G'-G) 



where the structure factor is 



SiG) 



iG 



(4) 



V^^ can be obtained in various way^^^EII, the em- 
pirical pseudopotential approach the atomic form factors 
are adjusted so the calculated energies at various high 
symmetry points fit experiment. In order for pseudopo- 
tentials to be transferable between polytypes (having dif- 
ferent G's)F^^(G) should be a continuous function of 
G. A wide vari ety of such model potentials have been 
used in literaturd^HMMI^ and use potentials of the form 



y^^(G) = (xiG + X2) [1 + exp(a;3G2 + Xi)] 



(5) 



where G = |G| x a/27r, and X j cLr6 adjustable parame- 
ters used to fit each material's cubic phase bandstructure 
to experiment. Model potentials that yield an accurate 
band structure of a known polytype should reliably pre- 
dict the band structure for the unknown polytype if the 
electronic environment in the two crystals structures are 
similar. 

We take spi n-orb it coupling into account by including 
the interactior{232ll 

{K',s'\Vso\K,s) = {K' xK)- {s'\rT\s) 



Xi p;(cos0k'.k)5(K'-K) 



(6) 

It is not necessary to expand Eq. [6] beyond / = 2 since 
group-IV semiconductors do not have core shells filled 
beyond d-orbitals. For Ge the terms up to ? = 2 in Eq. [7] 
are included, while for Si and diamond we only go up 
to Z = 1. Expanding Eq. |6]up to Z = 2, the spin-orbit 
coupling term is 



{K',s'\Vso\K,s) =^iiK' xK)-{s 
(Ap + ArfK'-K) 5(G'-G) 



(7) 



A; is a coefficient that can be written in terms of the core 
wave functions 



A, =Mi/3„;(K')/3„KK) 



(8) 



Pni{K)^C i'.jAn{2l + l)],,i{Kr)Rni{rydr (9) 
Jo 

where as are the Pauli matrices, K = G -|- k, 9 is the 
angle between K and K', and ^/ are adjustable pa- 
rameters used to fit the spin-orbit splitting energies to 
experiment^*. The overlap integral, /3„/, is constructed 
from the radial part of the core wave function, i?„i, which 
is an approximate solution to the Hartree-Fock equations 
which we obtain from Herman-Skillman tables^Sl. G is a 
normalization constant such that fini{K)/K approaches 
unity in the limit K goes to zero. Spin-orbit interactions 
are included for only the outer most p shells (n = 4 in 
Ge, n = 3 in Si, and n = 2 in diamond) and d shells 
(n = 3 in Ge). 



III. BAND STRUCTURES AND CRYSTAL 
SYMMETRIES 

A. Lonsdaleite Crystal Structure and Symmetries 

The LD /WZ crystal structure is constructed from two 
interpenetrating hexagonal close-packed (HCP) lattices, 
just as the diamond/ZB structure is constructed from 
two interpenetrating FCC lattices. For ideal crystals, 
the lonsdaleite lattice constant is related to the dia- 
mond lattice constant as a = acubic/V^, and the lat- 
tice constant along the c-axis (the [111] direction) is 
given by c = -\/l/u a. We assume an ideal LD crys- 
tal with u = 3/8 in this paper, giving a = 2.522 A 
and c — 4.119 A for diamond, both of which are in 
agreement with the experimental values of a = 2.52 A 
and c = 4.12 A?'. In the case of LD phase Si, a wide 
range of measured lattice constants have been repo rted 
with (a,c) = (3.84 A, 6.280 A)P, (3.84 A, 6.180 A)PS1, 
(3.837 A, 6.316 A^. Assuming an ideal LD crystal 
structure, we use a = 3.836 A and c = 6.264 A for Si. 
For Ge the ideal lattice constants are a = 3.993 A and 
c = 6.520 A, which is within 1 % of the experimental 
valuep2l of a = 3.96 A and c = 6.57 A. 

When viewed along the [111] direction, the inter-layer 
atomic bonds in LD lie in an eclipsed conformation, defin- 
ing the axis of hexagonal symmetry while in the CD 
structure the inter-layer atomic bonds are in a staggered 
conformation, making all four body diagonals of the cube 
equivalent. The nearest neighbors are the same in the 
two polytypes (due to their tetrahedral symmetry), and 
nine of the twelve next nearest neighbors are in identical 
positions in both crystals and the remaining three next 
nearest neighbors are equidistant. These structural simi- 
larities suggests that the local electronic environment will 
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be very similar in the two crystals, and hence the atomic 
form factors should be nearly identical in both polytypes. 

The LD crystal structure has a space group symme- 
try classification of _Dg^(or PG^/mmc). It has inversion 
symmetry in addition to all the symmetries of WZ. The 
irreducible representation of the space group of F are 
just the representations of the point group I?6/i (which 
has all the symmetries of Cgi, as well as inversion symme- 
try) . While moving along the kz direction the symmetry 
is lowered to Cqv The A, K and H all have point group 
symmetries of D^^^. D^h is isomorphic to C^y and is 
a symmorphic invariant subgroup of D^h ■ The L and M 
points have D2h symmetry. The point group operations 
must be followed by appropriate translations in order to 
obtain the irreducible representations of the wave func- 
tions at the high symmetry points. For example at the 
F point the point group operations must be followed by 
a translation t = (0, 0, c/2). 

The LD crystal structure has lower symmetry than the 
cubic diamond structure, and the SO interaction leads to 
additional lifting of orbital degeneracies. In the absence 
of spin-orbit coupling, the hexagonal crystal field of LD 
splits the p-like F15 state of cubic structure into a four- 
fold degenerate Fg and a doubly degenerate Fi. In terms 
of the p-orbitals, these states are pz — Fi are Px,Py — > 
Fg. With the inclusion of spin-orbit coupling, Fg„ splits 
into the Fg^ heavy-hole and the Fy^, light-hole. Therefore, 
all zone center states in LD belong to either F7, Fg, or 
Fg with either even or odd parity since LD has a center 
of inversion. Unlike WZ, there are no spin splittings in 
LD due to its inversion symmetry. 



B. Predicted Lonsdaleite Band Structures 

We use spherically symmetric local form factors fit to 
the CD band structure, which are then transferred to the 
LD crystal. The potentials should be transferable since 
the local electronic environments of the two polytypes 
are very similar. The lonsdaleite primitive unit cell has 
four atoms. We choose the origin, so that the atoms are 
located at vi = |ai + |a2, V2 = |ai + |a2 + ^as, V3 — 
vi + uaa and V4 = V3 -I- was. Where ai = (1, 0, 0)a, a.2 — 
(— 1, a/3, 0)a/2 and as = (0,0, c) are the the primitive 
lattice vectors. Substituting these atomic positions into 
Eq. [4j we obtain the following structure factor 



S 



1 



■ exp 



1 -f exp 



iGya 
iGxO, 



iGzUc 



GzUc 



(10) 



iGya 
2V3 



iGz 



where, Gj {j = x,y,z) are the components of the recip- 
rocal lattice vector G. 

The calculated band structure and the corresponding 
density of states(DOS) for diamond. Si and Ge in LD 
phase are shown in Figs. [l][3) The electronic band struc- 
tures are calculated in the irreducible wedge of the Bril- 
louin zone. The LD band structure is more complicated 




A R LUMZ FAA S HPK T F D(E) 

Figure 1. (a) Band structure for diamond in the lonsdaleite 
phase (b) Density of states. 
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Figure 2. (a) Band structure for Si in the lonsdaleite phase, 
(b) Density of states. 



its CD counter part due to its lower crystal symmetry. 
For a given energy range, there are roughly twice as many 
bands for the LD phase. The points, A and H, are spe- 
cial points where the energy levels stick together because 
the structure factor is zero there. 

Due to the similarities of the two crystals, many of the 
high symmetry points in the Brillouin zones of cubic and 
LD have a one-to-one correspondence with each other 
(just as in the case of ZB and WZ). This one to one 
correspondence is particularly useful in interpreting their 
respective band structures. 

The volume of LD's first Brillouin zone is about half of 
that of its cubic counter part. Therefore, if one were to 
take an intersection of the two Brillouin zones such that 
each of their F-points conincide, then the L-point in the 
cubic structure also coincides with the F-point in LB^. 



5 




A R LUMZ TAA S HPK T r D(E) 



Figure 3. (a) Band structure for Ge in the lonsdaleite phase, 
(b) Density of states. 



Thus, in the free electron model, the zone center LD/WZ 
energies can be directly predicted from the cubic T and L 
point energies. In single group notation (absence of SO 
interactions), the cubic Fi, Li and L3 states correspond 
to Fi, Fa and F5 respectively in LD. However, the pres- 
ence of the crystal potential will perturb the exact one to 
one correspondence of the high symmetry point energies. 
Because of this zone folding of the L- valley, indirect gap 
cubic materials with an L valley conduction band min- 
ima could be expected to have a direct gap in the LD 
phase (unless the energy of that state was significantly 
shifted by the crystal potential). 

There also are similar correspondences between the 
high symmetry directions of the two crystals. The 
A(F — >■ L) line in the cubic structure corresponds to the 
A(F A) line in the LEBU structure. Note that there 
are eight equivalent L directions in the cubic structure. 
Only the ones that are along the c — axis map on to the 
^hex point. The other six along with Xcub maps on to a 
point on the J/hex-line, two thirds away from the M^ex 
point. We label this point as M^^^. The Ac„b line maps 
onto a line joining Af^^^ and F/^ea;- The A^ub line is espe- 
cially important as the band gap for both diamond and 
Si lies along this line close to Xcub- 

We list the energies at various high symmetry points 
along with the corresponding irreducible representations 
of these states in tables [I] |III[ The irreducible representa- 
tions of the zone center states was determined by trans- 
forming the pseudo wave functions under the symmetry 
operations of the respective crystallographic point group. 
The zone center (F point) effective masses for diamond. 
Si and Ge for k parallel and perpendicular to the c-axis 
is shown in tabl dlVI 

From the band structure calculations, it is seen that 
in the LD phase, diamond and Si are indirect gap semi- 
conductors. Diamond has a band gap of about 4.767eV 



with band minima occuring at the ii'-valley. This is in 
agreement with earlier results obtained by Salehpour and 
SathpathjESl where the band minima was also shown to 
occur at K. Their estimated band gap from LDA cal- 
culations (with corrections) was about 4.5 eV. The LD 
band gap for diamond is also significantly smaller than 
its cubic phase band gap of 5.4 eV. 

In the case of the LD phase Si, the calculated band 
gap is 0.796 eV and the band minima is at the M- valley. 
The LD phase gap for Si is a lot smaller than its cubic 
badgap of about 1.1 eV. In our calculations, Gc is the 
only direct gap group-IV semiconductor in the LD phase. 
It has a Fg conduction band minima. The Fg symmetry 
is due to the fact that cubic phase Ge has a conduction 
band minima at the L valley, which folds over to Fg in 
the hexagonal phase. 

With the inclusion of SO interactions for LD, the top 
three valence states are typically (in order of decreasing 
energy) Fg, F7, F7 (normal ordering)''"' ^""^ or F7, Fg, F7 
(anomalous ordering) which results from a negative spin- 
orbit energy. In our calculations, the top three valance 
band states in Si and Ge have normal ordering in LD 
phase, while the top valance band states in diamond have 
neither normal nor anomalous ordering. Instead the top 
two valance band states for LD-diamond are F^ and F^ . 
The F^, r:^, states lie further bellow. 

The spin-orbit splitting energy. Ago and the crystal 
field splitting, Acr, can be extracted using the quasi cubic 
approximation which assumes the WZ/LD structur e to 
be equivalent to a [lll]-strained zincblende structur d^°^ l 
Aso and Acr are related to the F7„ light hole and F7„ 
light hole energies by 

i;(F^;^)-i;(Fg.)^- ^- + ^- 

±1(^{A,, + Acrr~^^y (11) 

where ^/u = a/c is 3/8 for an ideal WZ/LD structure 
assumed here. The band ordering and irreducible repre- 
sentations of the zone center states need to be identified 
before using this equation. We have complied a shorter 
table (tablejvj) for the LD phase semiconductors, that 
lists the band gap. Ago, Acr, and the offset between the 
valence band edges of each polytypes. 



IV. OPTICAL PROPERTIES 

A. Calculations 

At normal incidence, for any given polarization, 
the reflectivity acquires the simple form R — 
\{1 — rii) / {1 + ni)\^ , where, rii {i — x,y 01 z depending 
on the surface normal) is the complex index of refrac- 
tion. In the linear response regime, n{uj) = eiuj) and 
the complex dielectric function can be separated into real 
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Table I. Transition energies at various high symmetry points and the respective irreducible representations (IR) for the lonsdaleite 
phase of diamond. Note that the conduction band minima is in near vicinity of the K valley 
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r+ 


-1.7207 


Ar 


-4.3857 


is 


-7.6573 


Ms 


-6.321 


Hr 


-7.8319 


Kr 


-7.0031 


r+ 


-1.6793 


Ar 


-1.0552 


is 


-2.5751 


Ms 


-4.7847 


Hg 


-4.6241 


Kg 


-5.3748 


r+ 


-0.3496 


Ag 


-1.0322 


is 


-2.5516 


Ms 


-3.2709 


Hr 


-4.5975 


Ks 


-3.6509 




-0.0279 


As 


-1.0192 


is 


-2.5085 


Ms 


-3.01 


Hs 


-2.0413 


Kg 


-3.6304 







Ag 


-0.9962 


is 


-2.4906 


Ms 


-1.4971 


Hg 


-2.0303 


Kr 


-3.3739 


Ts 


1.4814 


As 


2.3662 


is 


1.2958 


Ms 


0.7957 


Hg 


2.3229 


Ks 


1.2191 




2.6962 


Ar 


2.3878 


is 


1.3036 


Ms 


2.9848 


Hs 


2.3319 


Kg 


5.7005 




2.7151 


As 


3.7479 


is 


2.4036 


Ms 


3.4078 


Hg 


4.1826 


Kg 


6.4501 


r+ 


3.4206 


Ag 


3.7614 


is 


2.4125 


Ms 


3.8669 


Hs 


4.1891 


Kg 


6.475 


r+ 
^ 7 


4.2135 


Ar 


3.7624 


is 


6.695 


Ms 


5.8173 


Hg 


6.8291 


Kr 


6.8974 


^ 9 


5.1671 


Ag 


3.7759 


is 


6.705 


Ms 


7.2281 


Hs 


6.8401 


Kr 


6.919 


r+ 

8 


5.1779 


As 


7.7044 


is 


9.1484 


Ms 


7.6228 


Hr 


7.3607 


Ks 


6.9399 



Tabic II. Transition energies at various high symmetry points and the respective irreducible representations(IR) for the lons- 
daleite phase of Si. Note that the conduction band minima is very close to the M valley 



IR 


E (eV) 


IR 


E (eV) 


IR 


E (eV) 


IR 


E (eV) 


IR 


E (eV) 


IR 


E (eV) 


r+ 
^ 7 


-12.4692 


Ar 


-11.7515 


is 


-9.9534 


Ms 


-10.5753 


Hr 


-9.5329 


Kr 


-9.7671 




-10.4783 


As 


-11.751 


is 


-9.9492 


Ms 


-8.9518 


Hs 


-9.4188 


Ks 


-9.5896 




-7.3386 


As 


-4.6049 


is 


-7.6576 


Ms 


-8.431 


Hs 


-8.0632 


Ks 


-7.471 


r+ 
^ 9 


-1.7267 


Ar 


-4.5941 


is 


-7.6503 


Ms 


-6.548 


Hr 


-7.8404 


Kr 


-7.3368 


r+ 


-1.4617 


Ag 


-1.1146 


is 


-2.4754 


Ms 


-4.3358 


Hr 


-4.1686 


Ks 


-4.787 


r+ 
^ 7 


-0.4896 


As 


-1.1069 


is 


-2.4673 


Ms 


-3.023 


Hg 


-4.0578 


Kg 


-3.3445 


r+ 
^ 7 


-0.1293 


Ar 


-0.8801 


is 


-2.3285 


Ms 


-2.6832 


Hg 


-1.9988 


Kr 


-3.3327 


Tg 





Ag 


-0.8723 


is 


-2.3199 


Ms 


-1.4302 


Hs 


-1.9369 


Kg 


-2.951 


r8 


0.3103 


As 


1.0951 


is 


0.7849 


Ms 


0.6828 


Hr 


1.7489 


Ks 


1.5267 


r+ 


0.7659 


Ar 


1.1023 


is 


0.7877 


Ms 


2.0997 


Hg 


1.8076 


Kr 


5.5903 


r| 


2.5368 


Ar 


3.538 


is 


2.4102 


Ms 


3.3933 


Hs 


4.1879 


Kg 


5.6603 




2.6714 


Ag 


3.5433 


is 


2.4151 


Ms 


3.7561 


Hr 


4.2033 


Kr 


5.6979 


r+ 
^ 7 


3.3102 


Ar 


3.6279 


is 


6.3148 


Ms 


5.0137 


Hr 


6.5626 


Kr 


5.7016 


r8 


4.7797 


Ag 


3.6332 


Lr, 


6.319!) 


Mr, 


6.6232 


H- 


6.6076 


Kr 


5.7153 




1.853() 


Ar 


0.08()2 


is 


8.;i.")ir) 


Mr, 


6.917 


H, 


6.G9G5 




G.1529 



Table III. Transition energies at various high symmetry points and the respective irreducible representations (IR) for the 
lonsdaleite phase of Ge. This is a direct gap semiconductor with a band gap of 0.310 eV. 
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diamond 


Si 


Ge 


m|| 


mx 


m|| 


mj_ 


m| 


m± 


1.1376 


1.1947 


1.1365 


1.1522 


1.1887 


1.1988 


0.1915 


1.1724 


0.2697 


1.1129 


0.3375 


1.1500 


0.1482 


1.8260 


0.184 


2.4867 


0.2153 


11.7707 


1.1656 


0.3104 


1.3479 


0.1843 


1.347 


0.1550 


1.1697 


0.3425 


1.3598 


0.1916 


1.4218 


0.1404 


0.1810 


0.6821 


0.1028 


0.7687 


0.0587 


0.3416 


0.3414 


0.2437 


0.5481 


0.201 


0.1484 


0.0871 


0.3418 


0.3226 


0.5637 


0.2128 


0.6035 


0.0672 


0.7860 


0.3785 


1.0483 


0.1224 


1.0563 


0.0852 


0.7877 


0.3588 


0.6496 


1.096 


0.0516 


0.0410 


0.2504 


0.9981 


0.6548 


0.9096 


0.6213 


1.0951 


0.2504 


2.6106 


0.1527 


1.0882 


0.6604 


1.5539 


0.9974 


0.5130 


4.0888 


0.1375 


4.5651 


0.4412 



Table IV. Zone center effective masses, parallel and perpen- 
dicular to the c axis, for diamond, Si and Ge. 





E, (eV) 


Aso (eV) 


A„. (eV) 


AEvB (cV) 


diamond 


4.7672 (Kg) 


0.0045 


1.6950 


0.0286 


Si 


0.7957 (Ms) 


0.044 


0.3336 


-0.1484 


Ge 


0.310 (Fg) 


0.404 


0.1730 


-0.1454 



Table V. Energies of the lonsdaleite phase of diamond. Si 
and Ge. The symmetry of the conduction band minimum 
is indicated in parenthesis along with the band gap. A^o 
and Acr are the spin-orbit spli ttin g and crystal-field splitting 



energies extracted using Eq. (hljl. AEvb = Eyb" - Eyg, 



is the energy difference between the top of the valance bands 
for the two polytypes. 



and imaginary parts e(cj) = + which are re- 

lated to each other by the Kramers-Kronig relations. 

All dielectric function calculations in this paper are 
carried out in the long wavelength limit, [i.e. assum- 
ing only direct band-to-band transitions (same k). We 
obtain t"{uj) using our empirical pseudopotential wave 
functions. In the electric dipole approximation, the di- 
rect transition between an initial state, I, and a final 



state, J, e"(w) is given by 



(12) 



BZ 



where is an integration over the entire Brillouin zone 
(BZ), is a sum over all initial valance band and final 
conduction band states, and ^^^(k) and i?c(k) are the 
valance and conduction band energies at their respective 
ks. The delta function is approximated by 

5{AE - nw) « 2 (1 -(- cosh[7(A£; - nuj)]y^ (13) 

where 7 is an adjustable damping parameter that can be 
used to phenomenologically incorporate lifetime broad- 
ening effects. We used 7 — 100 e V~^ which gives a tran- 
sition linewidth of about 35 meV'iSSI 



The real part of the dielectric function, e'{oj) is then 
obtained using the Kramers-Kronig relations 



duj' 



(14) 



where V is the Cauchy principle value. 

The momentum matrix elements, for band to band 
transitions, between the initial state / and and final state 
J is given by 



M/j(k) = {^imWjm) 



(15) 



where, p is the momentum operator. We calculate Mjj 
using the pseudopotential wave functions, which can be 
written as 



'/'/,k(i-) 



G 



c/(k,G)exp[i(k+G) -r] 



(16) 



where, c/(k, G) are the eigenvector coefficients, for the 
state, obtained by diagonalizing the pseudopotential 
Hamiltonian at a given wavevector k. Using Eq |16[ Mjj 
can therefore be rewritten in terms of the expansion co- 
efficients as follows. 



ij 



(k) =i^cl(k,G)cj(k,G)[(k-f G) 



(17) 



where, e is the polarization vector. The frequency de- 
pendent imaginary part of the dielectric function, e"(w), 
is then be calculated using EqjlTj for light polarized par- 
allel and perpendicular to the c-axis. In order to evaluate 
the Brillouin zone integral in EqfT2| we have used a set 
of 4.5 X 10"* special k points. These special k poin ts ar e 
generated using the scheme of Monkhorst and PadP^ 

The momentum matrix elements calculated using the 
pseudopotential wave functions need to be corrected for 
the missing core state^^. One way to do this is to include 
the commutator of the n onlocal pseudopotential and the 
position operatoj ^"^ * ^"^ *, while other pro pos ed methods 
involve the inclusion of a core repair ternPSI. However it 
should be noted that both techniques cause small changes 
to the dielectric function (typically less than 5%). Some 
first principles calculations have even shown that there is 
almost no difference between dielectric functions calcu- 
lated using ab initio pseudopotential wave functions and 
those calculated using true electron wave functionsSfiJ. 

Our efforts to calculate the dielectric functions of the 
LD phase of group-IV semiconductors however greatly 
aided by the fact that the cubic phase dielectric func- 
tions of diamond, Si and Ge are already well known. We 
take advantage of the fact that the pseudopotentials are 
being transferred between polytypes. Due to the simi- 
larity of the electronic environment around an atom in 
the two polytypes, corrections to account for the missing 
core states should be nearly the same and transferable 
between polytypes. This method also accounts for local 
field effects and static screening effects. 



First the dielectric functions for cubic diamond, Si and 
Ge are evaluated. We then normalize the calculated 
— 0) to the experimentally known static dielectric 
constant by making use of the optical sum rule 

2 f°° e"(uj)duj , , 

eo = 1 + K- / ^ ' 18 

where k is a scaling constant which is adjusted so that 
the calculated eo matches experimental results from Ref. 
mop . Next, the LD phase dielectric functions are cal- 
culated based on their respective band structures. The 
final LD phase e" and eJI are obtained by appropriately 

scaling their respective static dielectric functions and 

e'o- This scaled ej; and et! is obtained using the /—sum 



rule (Eq. 18 1 and the constant C (which is already de- 



termined from the cubic phase calculations). 

The static dielectric function needs to be corrected for 
the missing core states and other effects such as the in- 
fluence of local fields and static screening. The local 
field e ffects o ften arise from the inhomogeneity of the 
crystalliiJSI! and are difficult to account for in dielec- 
tric function calculations. However as they are known 
to alter the static dielectric function our scaling method 
includes the local field effects as well. This is further jus- 
tified by the fact that the similarity of the local electronic 
environment in both polytypes. 



B. Predicted Dielectric Functions 



As mentioned earlier, all zone center states belong 
to either Fy, Fg or Fg representations. The allowed 
interb and t ransitions for LD zone center states are as 
foUow^smH, For light polarized parallel to the c-axis 
(i?||), the optical selection rules only allow transitions 
between same symmetry states, i.e. Fy o Fy, Fg -O- Fg 
and Fg o Fg . For light polarized perpendicular to the c- 
axis(i5_L), the allowed transitions are Fy o Fy, Fg o Fg 
, Fg o Fy and Fg o Fg. Note, that the Fy o Fg transi- 
tion is forbidden for all polarizations. 

The other allowed high symmetry point transitions for 
E\\ are A,8 ^ ^7,8, Ag ^ Ag, ^4,5 o ^^4,5, Kq ^ Kq, 
lf4^5 o H^ c^ and iJg ^ Hq. Whereas for E± the allowed 
transitions are Ay^g -o- Ay^g, Aq o Ay^g, -^^4,5 ^ K^, and 
^4 5 ifg. All M and L valley transitions are allowed 
for all polarizations {i.e. M5 o M5 and L5 O L^^^^. 

In our calculations, Ge is the only group IV semicon- 
ductor that has a direct band gap in the LD phase. Its 
conduction band minima has Fg symmetry while the top 
of the valance band has Fg symmetry and hence LD phase 
Ge will only be optically active for E^. For E\\ only 
transitions from deep within the valance band will be al- 
lowed. The LD phase diamond and Si on the other hand 
are indirect gap semiconductors and will not be optically 
active. 

We calculate the frequency dependent e'(w) and e"(w) 
for parallel and perpendicularly polarized light for dia- 
mond. Si and Ge as shown in Fig. [4] As seen in Figj5] 
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hco (eV) ha (eV) 
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hco (eV) h(o (eV) 
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Figure 4. Real and imaginary parts of the complex dielectric 
function as a function of optical frequency shown here for 
light polarized {a)E± in lonsdaleite phase diamond, (b)i5|| 
in lonsdaleite phase diamond, (c)E± in lonsdaleite phase Si, 
(d) in lonsdaleite phase Si, {e)E± in lonsdaleite phase Ge 
and {f)E\\ in lonsdaleite phase Ge 

the corresponding reflectivity spectra for both polariza- 
tions shows several peaks which originate from interband 
transitions along various high symmetry points. Each of 
these crystals have distinct spectral features depend on 
the details of their electronic structure. For Ge and Si the 
most prominent features are typically seen up to about 
4 eV, whereas for diamond the most prominent peak is 
seen at about 12 eV. All the group-IV LD phase semi- 
conductors are optically anisotropic as expected. How- 
ever, diamond exhibits greater optical anisotropy than Si 
or Ge as seen in the reflectivity spectra (Fig[5|. Notice 
that around 12 eV, the absorption is significantly more 
for E± due the much smaller e± structures at the same 
frequency. 

In order to elucidate the the dielectric function's vari- 
ations about the fundamental absorption edge (FAE) , we 
fit the numerically calculated ej_ and e||, bellow the Re- 
strahlen band, to a classical Lorentz oscillator, whose real 
and imaginary parts are 




tico (eV) 




5 10 15 20 



tio) (eV) 




5 10 15 20 

hco (eV) 

Figure 5. Calculated reflectivity spectra at normal incidence 
for the lonsdaleite phase of (a) diamond, (b) Si and (c) Ge 



Material 


^0^ 


4' 


diamond 


5.7988 


6.3192 


Si 


12.1808 


11.8518 


Ge 


16.0066 


14.2896 



Table VII. Calculated static dielectric constant for light po- 
larized parallel and perpendicular to the c-axis, for the three 
lonsdaleite phase group IV semiconductors. 



Here, /, flj and Tj were used as fitting parameters and 
are listed in Table |VI[ Typically, / represents the oscilla- 
tor strength, Tj the relaxation rate and flj is a resonance 
frequency term. 

These fits in terms of the Lorentz oscillators provides 
an analytical expression for the dielectric function's dis- 
persion relations over the spectral region where there are 
no discontinuities. Therefore these dispersion relations 
are valid only within a certain cut-off frequency (ujk) 
as listed in Table |Vl]). The fitting parameters, such as 
/s, r^s and fi^s, depend on luk and do not necessarily 
represent trends in the optical properties of these semi- 
conductors. These analytic dispersion relations could be 
useful for modeling optical devices and multi-layer thin- 
film structures. 

The static dielectric constants were calculated 
from the imaginary part of the dielectric functions 
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(eo = 1 + 2 e"{uj)duj/TTUj), for light polarized parallel 
and perpendicular to the c axis, and are listed in Table 
|VII| In general, the semiconductor with the higher atomic 
number(Z) has the larger static dielectric constant. We 
also see that eQ < ejj in the case of diamond (which has 

a direct gap), whereas Eq > e[j for Si and Ge (which 
have Fg direct gaps). 

A simple explanation for this can be provided based 
on the optical selection rules. For a given material, Co 
depends on the number of allowed transitions and the os- 
cillator strength of each transition. The closer the bands 
are, the stronger the dipole transitions will be. Consider 
a small region of u around the direct band gap where the 
transitions are the strongest for zone center states. For 
E± a transition between the Fg heavy-hole (HH) and 
the Tj (or Fg) conduction band is allowed. All three 
materials Si, Ge and diamond will be optically active. 

f 1 2) 

Whereas for the Fy ' ' light-hole/split-off hole to Fy 
conduction band is only allowed for diamond. Si and 
Ge are completely optically dark in this case. Hence for 
diamond has more allowed transitions than for E± 
(assuming similar oscillator strengths) and therefore its 

Sq < ejj . Whereas the opposite is true for Si and Ge. 



V. SUMMARY 

We have calculated the electronic band structures and 
dielectric functions for diamond, Si and Ge in lonsdaleite 
phase using transferable model potentials, including spin- 
orbit coupling. The potentials should be accurate since 
the local electronic environment for the cubic and lons- 
daleite polytypes are very similar. It is seen that while 
diamond and Si remain indirect in the lonsdaleite phase, 
Ge is transformed into a narrow direct-gap semiconduc- 
tor due to zone folding effects. Hence LD-type Ge will 
be optically active, which could make it extremely useful 
for technological applications. We have also tabulated 
a number of parameters such as high symmetry point 
energies, their irreducible representations and effective 
masses, which could be useful for constructing k ■ p type 
models. We calculated the frequency dependent complex 
dielectric functions up to 20 eV for light polarized paral- 
lel and perpendicular to the c-axis in the dipole approx- 
imation. We find strong optical anisotropy, making LD 
phase materials potentially useful as nonlinear crystals 
since their optical birefringence enables them to satisfy 
phase matching conditions. 



VI. APPENDIX 



In this appendix we provide the details for our band- 
structure calculations for the cubic phase of diamond. Si 
and Ge. The form factors of the cubic polytypes are re- 
quired for obtaining the bandstructures of the hexagonal 




Figure 6. Local atomic form factors(V^) for diamond(C), Si 
and Ge. Note that the y-axis for diamond is on the right. 



diamond 



10 

% 
1-10 

LU 

-20 




u;Ki r 



Figure 7. Calculated band structures, with the inclusion of 
spin-orbit interactions, for the cubic phase of (a) diamond (b) 
Si and (c) Ge 



polytypes. The Hamiltonian, given in Eq T] is diagonal- 
ized with a plane wave basis cutoff of |G < GAir/a. As 
a starting point, our continuous form factors for Si and 
Ge were first fit to the discrete form factors of Ref.[98] 
and for diamond, they were fit to the discrete form fac- 
tors of Ref. |114| . The adjustable parameters Xi, and 
in Eqs. [5] and [8 were then adjusted to fit the calculated 
band structure to the experimental energies of the band 
extrema of the cubic materials. A modified simulated 
annealing method was used for the fitting procedure (see 
RefE for more details). 

Five adjustable parameters for diamond and Si, and 
a total of six adjustable parameters for Ge were used 
to fit the calculated bandstructure to seven different ex- 
perimental transition energies (which were obtained from 
Ref. I110|) . Additional constraints were imposed to ensure 
the correct band ordering of valence states and conduc- 
tion states. 

As can be seen in tablejVinj our fits are in excellent 
agreement with experiment. The fitting parameters for 
the form factors, used for obtaining the results shown 
in table I Villi are given in table |IX| and the atomic form 
factors themselves are shown in fig|6] The calculated 
ZB bandstructures of these group-IV semiconductors are 
shown in fig(7] Note that for cubic diamond and Si, the p- 
like conduction band states are bellow the s-like state. In 
addition, these are also indirect gap semiconductors with 
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diamond 


Si 


Ge 


Transiton 


Expt. (eV) 


Calc. {eV) 


Transition 


Expt. (eV) 


Calc. {eV) 


Transition 


_Ba;pf. {eV) 


Calc. (eV) 


'-7c 8d 


7.3 


7.3 


7c 8i> 


3.35 


3.35 


v~ — r+ 

'- 7c 8u 


0.898 


0.898 


-"-Be 8d 


12.9 


12.854 


'-lie Sv 


4.15 


4.148 


'-6c '- 


3.22 


3.223 


r+ -r+ 

'- Sv '- 6v 


26 


26.779 


8u -"-61) 


12.5 


12.707 


8v -"-611 


12.6 


12.286 


^6c ^(,v 


10.5 


10.508 




2.05 


2.072 


^6c 8 


0.76 


0.76 


X^c — X^v 


12.9 


12.89 


X^c — Tg 


1.13 


1.13 




1.16 


1.081 


Aso 


0.006 


0.006 




0.441 


0.441 




0.297 


0.297 


Ko 






Ko 








0.2 


0.2 



Table VIII. A comparison between targeted experimental transition energies taken from Ref.jTTOj and the calculated values for 
cubic diamond, Si and Ge. For each material, the first column shows the targeted transitions, the second column shows the 
targeted energies and the third column shows the converged results from fitting the pseudopotentials. The spin-orbit splitting 
energies are A^o = E^„ — E^^ and for Ge, A!,^ — — 



Material 


XI 


X2 


X3 


Xi 


IJ-l M2 


diamond 


444.305 


-1716.53 


0.0263 


6.2294 


12.515 


Si 


8.2808 


-54.1842 


0.0116 


4.7922 


0.0536 


Ge 


0.0791 


-0.5737 


0.0247 


-1.2269 


0.2413 8.4571e-17 



Table IX. The XjS are the fitting parameters for the form factors and spin-orbit splitting energies and /ii and fj,2 are the fitting 
parameters for the spin-orbit coupling. Note that the form factors are in units of Ry. 



the conduction band minima lying in very close proximity to the X-valley along the A-direction. 
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